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VjQ , Single top-quark production, via weak-interaction processes, is an im- 

00 ■ portant test of the standard model, potentially sensitive to new physics. 

However, this measurement is much more challenging at the Tevatron than 
nI originally expected. We reexamine this process and suggest new methods, 

using shape variables, that can supplement the methods that have been 
discussed previously. 



The electroweak production of single top quarks is an important stan- 
dard model process which the Tevatron is guaranteed to be able to study. 
This reaction, which has been investigated previously Q][2|, is interesting 
<-j ! both because it provides a direct measurement of the Vtb CKM element 

and because it is sensitive to deviations of top quark physics from standard 
model predictions [3]. Limits on single top production from Run I at the 
Tevatron have been published [3J |^j, and the first Run II limits have ap- 
peared [HI Ej. This discussion will provide a brief overview on our recent 
work on this subject [S]. The executive summary is that a) this process will 
be more difficult to study experimentally than previously thought, b) we 
have developed new analysis methods, which will improve the significance 
of the measurement, and c) even with our new analysis methods the un- 
ambiguous detection of this process remains a considerable challenge. In 
particular, we suggest making use of more of the information encoded in 
the shape of the signal, in a way that will be less sensitive to systematic 
errors than a simple counting experiment. On the other hand, we also show 
that the size of and the uncertainties in the VU-plus-jets background cause 
serious problems that at present make the measurement difficult at best. 

At Fermilab energies, the u tb" production of a top quark and bottom 
antiquark by an s-channel W boson, as can occur through the diagram in 
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Fig. 1. Single top quark production via a) an s-channel W, b) a t-channel W. 



Fig. El), has a lower cross-section than "tbq" production via a t-channel W 
boson of a t, b, and an extra quark jet near the beam axis, as occurs through 
diagrams such as that in Fig. 

From Fig. ^>) one can see that a tbq event has a final partonic state 
consisting of at least the following: a charged lepton, a neutrino, b and 
b quarks, and a light quark. Thus, in selecting t-channel signal events, 
one asks for (a) one 6-tagged jet, (b) significant missing energy (from the 
neutrino), (c) one and only one isolated or /i ± , (d) at least one non-6- 
tagged jet. Typically, the highest-pr non-6-tagged jet in a t-channel event 
is that from the light quark. Also, in typical events a t quark can be 
reconstructed from the tagged jet and a W, itself reconstructed from the 
charged lepton and the missing energy. The tb process has a b quark jet 
and a b jet, along with a lepton and a neutrino. In a significant number of 
events, one of the two 6-jets will not be tagged, so that the same criteria 
used for tbq will have moderately high efficiency for this process as well. 

The main backgrounds to single top production, which all can imitate 
the signature just described, are (1) "it", top quark pair production, pri- 
marily from events where only one of the top quarks decays leptonically; (2) 
"QCD" events with fake electrons or muons, or with real muons from heavy 
flavor decays; and, most problematic, (3) a Wj n " events with a leptonically- 
decaying W plus some number n > 2 of quark or gluon jets. Processes of all 
three types occur at rates well above the signal rate and so we must provide 
specialized cuts and other methods to reduce their contribution to the final 
event sample. For example, the ti background events tend to exhibit consid- 
erably more transverse energy than true single top events, to have more jets, 
and to be more spherical. Experience with the QCD background suggests 
that it is generally reduced to a level comparable with the signal [7j by fairly 
standard cuts and we will not focus on it here, except to note that our event 
shape methods will tend to also substantially reduce the QCD component of 
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the event sample. The VK-plus-jets background is much more complicated. 
The Wj n events potentially entering the sample consist of a real W boson 
decaying leptonically, and at least two other quarks or gluons in the final 
state. While the Wj n events do tend to have smaller energetics than single 
top, and, of course, have no reconstructible t quark; the number of events is 
so large, and the energy resolution at Fermilab is sufficiently broad, that the 
Wj n events form a large and problematic background to single-top produc- 
tion. Reducing the systematic error on the prediction and/or measurement 
of this process is essential for success. 

To model both signal and background, we have used MadEvent [Hj to 
generate events, Pythia jlUj to then simulate showering and hadronization, 
and PGS 11 to act as a fast detector simulation. (See 8 for the details of 
the calculations.) Jets containing a b quark (either perturbatively or pro- 
duced during showering) are taken to be tagged with an efficiency of the 
form 0.5 tanh (py/36 GeV), where px is the transverse momentum of the jet. 
Jets containing a c quark (either perturbatively or produced during shower- 
ing) are taken to be tagged with a rate of the form 0.15 tanh (pt/42 GeV); 
while jets containing no heavy flavor are taken to be mistagged with a rate 
of the form 0.01 tanh (p T /80 GeV). 

We have chosen for this study to cut also on the quantity 

H T = ^J(pt)i + {vt)i +$t , 

jets 

which, compared to the signal, is smaller for Wj n and larger for ti. In this 
expression the sum is over all jets with px > 20 GeV and \r]\ < 3.5, (pr)i 
is the magnitude of the transverse momentum of the i th jet, {pt)i that of 
the lepton, and-^r is the missing transverse energy in the event. Since the 
signal involves a t-quark, we also impose a requirement that the invariant 
mass of the lepton, neutrino, and the leading tagged jet be within a window 
around the top quark mass. 

We find that such cuts cannot decrease the backgrounds to the point 
that they are comparable to the signal. Two choices of "intermediate" and 
"hard" cuts are indicated in Table ^ The resulting numbers of expected 
events for an integrated luminosity of 3 fb _1 , summing over and (and 
thus including both t and t), can be seen in Table El Consistent with pQ, we 
find that while all of the cuts contribute to the background reduction, the 
Wj n channel is reduced primarily by a combination of the stiffer p? cuts 
and the "mt" cut, while the ti background is affected primarily by the upper 
Ht cut. While the intermediate cuts reveal a signal to background ratio of 
1:7.4, this improves to 1:4.9 using the hard cuts. (The basic cuts used as 
a starting point in [H] yield a signal to background ratio of approximately 
1:21.) This result is, at best, disappointing and more pessimistic than 
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Item 


PT 


V 


PT 


V 


e ± 


> 15 GeV 


< 2 


> 15 GeV 


< 2 


MET 0) 


> 15 GeV 




> 15 GeV 




jet (6-tag) 


> 20 GeV 


< 2 


> 60 GeV 


< 2 


jet (no 6-tag) 


> 20 GeV 


< 3.5 


> 30 GeV 


< 3.5 




Min 


Max 


Min 


Max 


Hj 1 


> 180 GeV 


< 250 GeV 


> 180 GeV 


< 250 GeV 


"m t " 


> 160 GeV 


< 190 GeV 


> 160 GeV 


< 190 GeV 



Table 1. Representative "intermediate" cuts (columns 2 and 3) and "hard" cuts 
(columns 4 and 5). 



Channel 


"Intermediate " Cuts 


"Hard" Cuts 


tbq 


67 


30 


tb 


27 


13 


tt 


140 


57 


Wjj 


550 


152 


(tbq + tb) 1 (tt + Wjj) 


0.14 


0.21 



Table 2. Numbers of events for 3 fb 1 (summed over t and t, e and fi channels) for 
the two sets of cuts in Tabled 

(For a detailed comparison with the results of see |S].) We are concerned 
that systematic errors in the understanding of the background will plague a 
direct counting experiment at a level that will make any claims of discovery 
suspect. 

Under the assumption that a counting experiment is indeed insufficient, 
we turn to observables that (as in [31 El ) make use of other aspects of the 
signal. For the dominant tbq production process the light quark, forward 
jet is a distinctive signature, which the backgrounds do not share. We make 
use of this feature in our analysis of both signal and background events by 
including in our event definition the properties of "the jet" , defined to be the 
highest-pT non-6-tagged jet. Since the tbq process arises from an initial state 
light quark or antiquark scattering off a gluon, the tbq system is typically 
boosted in the direction of the initial quark. Moreover the motion of the 
final state light quark is typically in the same direction as its parent quark 
(and the proton). (The reverse is true for t production.) The structure of 
the electroweak interactions in both the production and decay of the top 
quark tends to align the momentum direction of "the jet" with that of the 



ISMDsde printed on February 2, 2008 



■5 



final-state charged lepton. Thus the momentum vectors of the lepton and 
"the jet" are correlated as a result of both kinematics and spin polarization 
effects. 

The tt process shares none of these properties. At tree-level, tt is sepa- 
rately C and P invariant, and shows none of the above asymmetries. The 
parity-invariance is violated at the next order ^2], at a level both small 
and calculable. Moreover, there should be no strong correlation between 
the momenta of the lepton and the jets in the final state. Similar consider- 
ations apply, to a good approximation, to those QCD events which might 
pass our cuts, mainly events with fake leptons or with isolated leptons from 
heavy flavor decays. (See jS] for additional discussion of this point.) 

The asymmetries and correlations in Wj n are similar to those of tbq 
(unfortunately) although less pronounced. As with a t quark, a W + is 
most likely to be produced moving in the proton direction. This leads to 
parity asymmetries which, though relatively small, are still quite large in 
absolute size compared to the signal. Correlations between the lepton and 
"the jet" are nonzero, though relatively small. Unfortunately, the size of 
the asymmetries and correlations in Wj n appears to be very sensitive to 
assumptions, cuts, Monte Carlo parameters, and tagging, and will be a 
source of significant systematic error. 

To make the use of these special properties of the signals and back- 
grounds we consider two-dimensional distributions in pseudo-rapidity. Since 
the pp initial state of the Tevatron is a CP eigenstate, the differential cross- 
section d 2 a + /dr]jdr]£ describing the rapidity distributions of "the jet" and 
the charged lepton in events with a positively charged lepton, and the corre- 
sponding distribution d 2 a~ /drjjdrje for processes with a negatively charged 
lepton, must be CP invariant: 

d 2 a + , . d 2 a~ . 
dr]jdr]e drjjdrji 

Consequently, we can combine data from positively and negatively charged 
leptons by defining a lepton- charge-weighted pseudo-rapidity, f)j = Qirjj, 
fje = Qir]i, where Qi is the lepton charge. The remainder of this discussion 
is based on the explicitly CP-invariant differential cross section 

To illustrate the characteristic properties of the various channels we will 
study a simulated event sample defined by cuts somewhat more "relaxed" 
than those of Table Q keeping a larger fraction of both the signal and the 
background. The simulated differences in shape are summarized in the 
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Fig. 2. Differential cross-section (cPa/dfijdfji, summed over t and t, e and /i) for 
the a) tb channel b) tbq channel, c) it channel, and d) Wjj channel. 



contour plots of Fig. EJ which give the distributions (d 2 a / 'dfjjdfji) of the 
various processes in the (fjj,fje) plane. 

Figures^ and[2t illustrate the relative symmetry of the tb and it chan- 
nel, to be compared to the small asymmetry in the Wjj channel in Fig. [2ji 
and the large asymmetry in the tbq signal channel of Fig. |2j3. In particular, 
the tbq signal lies dominantly in the B quadrant, where the four quadrants 
of the [fjj,fji) plane are labeled in the form ^ | ^ . Meanwhile tb is smaller 

with a less distinctive shape. By contrast, ti is large though symmetrically 
distributed, and Wj n is large but has small asymmetries between the four 
quadrants, with slightly more cross-section in the A and B quadrants than 
in C and D. 



To capture quantitatively these differences in shape between signal and 
background, we define three orthogonal functions in the {fjj,fji) plane. Any 
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differential cross-section can be written as a sum of three components 
d 2 a 

(fjj , fje) = F(rjj ,r)t)+F+ (fjj ,fj e ) + F- (fjj , fje) , 



dfjjdfje 



(1) 



where the components are of the form 

1 



F (fjj, fje) = 



d 2 a d 2 a » , 

(Vj,Ve) + -77. {~Vj, -m) 



dfjjdfje 



dfjjdfje 



d 2 a • d 2 a , „ 

+ ,~ w< + » ,~ i-Vj,Ve) 



dfjjdfje 



dfjjdfj e 



(2) 



d 2 



a 



dfjjdfje 



dfjjdfje 



dfjjdfje 
d 2 a 

dfjjdfje Vj,7]e 



Vj,-Ve) 



(3) 



F_(fjj,fje) = 



d 2 



a 



dfjjdfje 



d 2 



a 



dfjjdfje 



7-(-fjj,-Ve) 



(4) 



These functions are orthogonal in the sense that the integral of any (non- 
identical) pair over any symmetrically-defined region of the (fjj, fje) plane 
vanishes. Most importantly, the functions provide important physical in- 
formation about the symmetry properties of the differential cross-section. 
F and F + are parity-even while F_ is parity-odd; thus F_ = (within 
statistics) for any parity-even distribution. Meanwhile F + will also vanish 
if the distribution is parity-even and the leptons and jets are uncorrelated. 
By construction F and F + have four-way symmetry under reflection in the 
(fjj, fje) plane, while F_ has two-way symmetry (quadrants A and D are 
related, as are B and C, but quadrants A and B are independent). 

The fact that the signal has strong asymmetries and correlations not 
present in the backgrounds yields the following expectations, 



ptt > F tt . 



F a ~ p Wjn > F tbq , 



> F^ 



. ptbq 



F™ ~ 



+ 



ptbq 



F_ < F_ 



Wj n 



ptbq 



ptbq 



Thus the backgrounds are very large only in F, while the signal has a much 
larger role to play in the other functions, especially away from the center of 
the (fjj, fje) plane. 
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Fig. 3. Contour plots for the function F (fjj, fjt) for a) tb, b) tbq, c) tt, and d) Wjj 
channels (summed over t and t, e and fi) . 

These expectations are illustrated in Figs. 1210 and |SJ where the functions 
F,F + ,F_ are shown, for tb, tbq, ti, and Wj n . Note, in particular, the 
presence of the symmetry properties discussed above. 

Thus we are led (see [S] for more details) to suggest the following gen- 
eral approach. First construct, for both the data and the Monte Carlo 
simulations, the F, F + and F_ functions. Using these functions, as well as 
information obtained from other measurements, we systematically test our 
understanding of each process. The F function allows a measurement of 
the sum of the backgrounds without much contamination from signal. As- 
suming that the separation of Wj n from tt can be obtained using the fact 
that the ti process can be measured and predicted with reasonable accuracy, 
we can then cross-check our understanding of the Wj n background using 
part of the F- distribution (located roughly in quadrant A) where the signal 
is negligible. Finally, one can attempt to measure the signal from F + and 
from another part (located largely in quadrant B) of the F_ distribution. 

A complete analysis of the above ideas has been carried out in jS], in- 
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Fig. 4. Contour plots for the function F + (f)j,fj(,) for a) tb, b) tbq, c) tt, and d) Wjj 
channels (summed over t and t, e and fi) . 



eluding estimates of the associated uncertainties, both statistical and sys- 
tematic. The primary lesson is that, even when using the shape variables 
discussed above, the contributions of the backgrounds, especially Wj n , can- 
not be rendered truly small compared to the signal. Thus it is essential that 
we understand the backgrounds with high precision, including their shapes 
as employed here. This last point is especially challenging for the Wj n case, 
which receives comparable contributions from events with tags of b quark 
jets, tags of c quark jets and mistags of jets with no heavy flavor. Like- 
wise there are comparable contributions from when the heavy-flavor quark 
is present in the short-distance scattering and when it is produced in the 
subsequent showering. This means that our ability to accurately simulate 
the background depends on our understanding of a multitude of different 
specific channels, each with their own shapes and their own dependence on 
cuts, parton distributions and other parameters. 

In summary, a counting experiment for discovery of electroweak single- 
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Fig. 5. Contour plots for the function F_ {f)j,fji) for a) tb, b) tbq, c) tt, and d) Wjj 
channels (summed over t and t, e and n) . 



top production appears very challenging. We have explored the possibility 
of using the distinctive shape of this process to separate it from background. 
We have noted that the distributions for tt and for QCD in the (fjj,fj() plane 
are largely symmetric, while that of the signal is not; the Wj n (W-phis- 
jets) background is intermediate between them. Constructing the functions 
F,F±, which have various symmetry properties, we find that the statistical 
and systematic errors in the functions F± , which are orthogonal to the func- 
tion F (that would be used in a counting experiment), can be brought close 
to reasonable size without using extreme cuts. Our method largely removes 
tt and QCD events from the observables, making extreme cuts on tt unnec- 
essary, and focusing attention on Wj n as the main background. The most 
challenging problem is understanding the shape of the Wj n background. 
We believe this calls for a dedicated study of the the rates, shapes, and fla- 
vor content (especially of bottom versus charm) of both Wj n and Zj n , with 
zero, one and two tagged jets, blended with theoretically precise predictions, 
and careful tuning and cross-checking of Monte Carlo simulations. 
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